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Abstract 

We describe the influence of electromagnetism on the phenomenology of if — > 
7T7T decays. This is required because the present data were analyzed without 
inclusion of electromagnetic radiative corrections, and hence contain several 
ambiguities and uncertainties which we describe in detail. Our presentation 
includes a full description of the infrared effects needed for a new experimental 
analysis. It also describes the general treatment of final state interaction 
phases, needed because Watson's theorem is no longer valid in the presence 
of electromagnetism. The phase of the isospin-two amplitude A2 may be 
modified by 50 — > 100 %. We provide a tentative analysis using present data 
in order to illustrate the sensitivity to electromagnetic effects, and also discuss 
how the standard treatment of e'/e is modified. 



I. INTRODUCTION 



In this paper we address the effect of electromagnetism on the phenomenology of K — > 
7T7T decays. Our previous work of Refs. has dealt mainly with the determination of 

structure dependent EM effects on the K — > titt amplitudes. It is the aim of the present 
work to attempt a complete phenomenological analysis. We start by briefly reviewing the 
standard treatment of K — > tttt decay amplitudes, pointing out that there is room for 
potentially important isospin breaking effects. We then focus on the effect of electromagnetic 
interactions (EM); we enumerate the main new features due to EM and describe their impact 
on the parameterization of the K — > tttt decay amplitudes. Our quantitative analysis begins 
in Section |J where we take up the problem of the infrared divergences which are induced 
by electromagnetism. We provide a complete description suitable for use in an experimental 
analysis in Section |T|. We also point out that the data on K —>■ titt do not lead to a direct 
extraction of the strong phase shift difference 5q — 62, because electromagnetism changes the 
rescattering phases of the amplitudes. Already a perturbative calculation has provided clues 
for the presence of sizeable effects ||. To account for this effect more generally, we provide in 
Section [TV] a suitable extension of Watson's theorem, obtained after writing and solving the 
unitarity constraints in presence of isospin-breaking interactions. Our goal throughout the 
paper is to relate the theoretical and experimental issues of these decays, in the hope that 
a future experimental analysis will be undertaken to resolve the substantial experimental 
uncertainties. 

The reason why present experimental information is not adequate is that most of the 
data was analyzed without the inclusion of electromagnetic radiative corrections. This 
means that the data in the Particle Data Tables is not fully reliable. Moreover, for certain 
quantities (namely 5q — 62 and the 1 = 2 amplitude A 2 ) this uncertainty is enhanced by 
a AI = 1/2 enhancement factor of 22. In Section [V] we demonstrate these effects by 
giving an illustrative data analysis, trying our best to interpret the present data. This 
step is necessarily tentative and likely partially incorrect, as it requires knowledge of the 
experimental procedure used to deal with soft photons when measuring the branching ratios. 
In the absence of detailed information from the Particle Data Group (PDG), we adopt the 
simplest theoretical framework, not necessarily corresponding to the real experimental setup. 
Within this simple framework we are able to show that the extraction of EM free quantities 
is quite sensitive to the treatment of infrared photons. 

The other interesting phenomenological issue concerns the effect of electromagnetism on 
CP-violating observables. In Section [VT| we discuss the impact of our findings on theoretical 
predictions for e'/e. In particular, we provide a quantitative estimate for the parameter 
f2 EM , the effect of a AI = 5/2 interaction, and the impact of the new rescattering effects on 
the phase of e'. We conclude the paper with a summary of our findings in Section VII. 



1 



A. Standard Phenomenology for K — > tttt Decays 

Let us start from the conventional phenomenological analysis of the decay amplitudes. 
There are three physical K — > nn decay amplitudes^] 

AkO^tt+tt- = A + - , A K °-*ir°K° = .4 00 , Ak+^tt+tt = -4+o ■ (1) 

We consider first these amplitudes in the limit of exact isospin symmetry and then identify 
which modifications must occur in the presence of electromagnetism. In the 7 = 0,2 two-pion 
isospin basis, the physical amplitudes are parameterized as: 

A + - = A e lSo + ] j^A 2 e^ , 

Aoo = A e i5 °-V2A 2 e iS * , (2) 
A +0 = ~A 2 e iS > . 

In the above, Aq^ are the AI = 1/2, 3/2 transition amplitudes corresponding to the tttt final 
states with isospin equal to 0, 2 respectively. They are real in the limit of CP conservation. 
The Si are the I = 0, 2 tttt scattering S-wave phase shifts at center of mass energy equal 
to the kaon mass. They enter the parameterization as prescribed by unitarity (the Fermi- 
Watson theorem). Knowledge of the experimental branching ratios [|J allows one to use 
Eqs. (0) to extract the isospin amplitudes. Neglecting the small CP- violation effect, we 
findQ 

A = (5.458 ± 0.012) x 10" 7 M^o , 

A 2 = (0.2454 ± 0.0010) x lO~ 7 M K o , (3) 
5 - 5 2 = (56.7 ±3.8)° . 

Throughout we express the K — > tttc amplitudes in units of 10~ 7 M^o, with M K o = 0.497672 
GeV. 

A careful inspection of these phenomenological results reveals some inconsistencies with 
other existing pieces of phenomenology and theoretical analysis. These clues seem to suggest 
that isospin breaking effects (like EM) can play an important role in the phenomenology 
of K —>■ tttt decays. The considerations we present below apply to all isospin breaking 
interactions. In this category one includes both EM and strong isospin breaking, produced 
by the difference in the up and down quark masses. In this work we are concerned exclusively 
with EM effects (also analyzed in Refs. [^-0]). For treatments of strong isospin breaking 
effects see Refs. |||§. 



*The invariant amplitude A is defined via ou t (7T^|^) in = i(27T) 4 S^( Pont - Pin ) (iA). 

^Knowledge of the phase difference 5o — 5 2 from other determinations poses the constraint cos(<5o — 
^2) > 0, implying that AqA 2 > 0. In this paper, we take Aq and A 2 as positive numbers. 
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Isospin breaking interactions will in general mix the amplitudes Aq and A 2 , thus generat- 
ing potentially big corrections to A 2 , proportional to A ■ a /it. A related issue concerns the 
presence of a AI = 5/2 component in the interaction. This problem has recently received 
attention in Ref. fpT0| . A AI = 5/2 component will distinguish between the amplitudes A 2 
entering in the K° and K + decays. The expression of Aq,A 2 ,A 2 in terms of Aaj is given 
by: 

^0 = A 1/2 , 

A 2 = A 3/2 + A b/2 , (4) 
At = A 3/2 - 2/3 A 5/2 . 

We expect the dominant AI = 5/2 effect to arise by combining the large AI =1/2 weak 
interaction with the AI = 2 component of the electromagnetic interaction. A combination 
of the AI = 3/2 hamiltonian with the AI = 1 interaction proportional to m u — is also 
expected to contribute. However, its effect is expected to be doubly suppressed (by the 
AI =1/2 rule and the smallness of m u — rrid)- 

A further problem with the isospin analysis is revealed by looking at the extracted 
phase shifts. The value 5q — 5 2 = (56.7 ± 3.8)° obtained from kaon decay data has to be 
compared with information coming from other sectors of low energy phenomenology. In 
particular, the value extracted from a dispersive treatment of im scattering data is So — S 2 = 
(45.2 ± 1.3±f;l)° fLl, and the prediction of ChPT is 5 - S 2 = (45 ± 6)°. These two 



determinations are mutually compatible. However, there is a sizeable discrepancy with the 
result obtained in the fit to K — > tttt. This can be ascribed to isospin breaking and to a 
non- vanishing A 5 / 2 . 

The above mentioned issues also affect a proper theoretical understanding of the direct 
CP- violation parameter e'/e. In particular, the leakage of the octet amplitude in A 2 (due to 
isospin breaking effects) brings an extra contribution to the CP- violating phase of A 2 . In 
the literature only the leakage due to m u — isospin breaking effects has been analyzed, 
and is found to be numerically important |8||| . Moreover, the presence of a AI = 5/2 
amplitude introduces an extra term in the usual formulae for e'. Finally, understanding the 
issue regarding the phase shift 5q — S 2 will provide a better theoretical determination of the 
phase of e'. 

The above considerations call for a careful analysis of electromagnetic effects on K — > tttt 
decays. 



B. Electromagnetism and the K — > -kit Amplitudes 
One can summarize the effects of electromagnetism on K — > tttt amplitudes as follows: 

1. First of all one has to deal with universal infrared (IR) effects, due to photons of long 
wavelength. These effects are common to all processes with charged external particles 
and do not depend on details concerning the original interaction. They are represented 
diagrammatically in Fig. [I], where the dark blob is seen as pointlike by the infrared 
photons. This class of contributions provides a Coulomb final state interaction (FSI) 
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(a) (6) (c) 



FIG. 1. Electromagnetic contributions relevant for the infrared behavior. 

phase and gives rise to IR divergent amplitudes. Such infrared divergences have to 
be canceled by considering the effect of soft real photons (Fig. |T](b)). As we shall 
show in Sections || and |V|, these effects can be described as modifying the phase 
space factor rather than producing effects on the amplitudes themselves. In order to 
perform an accurate phenomenological study, it is important to include this class of 
radiative corrections in the experimental analysis of the branching ratios (See Sect. ^ 
for details). 

2. There are structure dependent effects, sensitive to the form of the original interaction. 
These are hidden in Fig. |I] within the large dark vertices. We consider only corrections 
induced by the dominant (octet) part of the weak hamiltonian. These produce shifts 
in the isospin amplitudes, and are responsible for possible large contaminations of Ai. 
They also generate a AI = 5/2 amplitude. We have calculated these effects in ChPT 
and in dispersive matching. 

3. Finally EM affects the final state interaction. As a consequence the unitarity relations, 
determining the rescattering phases, are altered. The main modifications are due to 
the opening of the tttt+wj intermediate channels and the possibility of mixing between 
two-pion states in isospin 1 = and 1 = 2. These new features imply modifications 
to the unitarity parameterization, governed by an extension of Watson's theorem that 
we shall discuss at length in Sect. IV. 

Clearly the items enumerated above are intertwined and affect in various respects the 
way one analyzes the K — > tttt amplitudes. It is important to note that - although being per- 
turbatively small - the new interaction considered breaks the original isospin symmetry on 
which the parameterization of K — > tttt amplitudes is based. Therefore, to perform a com- 
plete analysis of EM effects to K — > tttt decays, one must understand how to parameterize 
the above mentioned effects. 



II. INFRARED BEHAVIOR AND ISOSPIN AMPLITUDES 

A. Defining Infrared Finite Amplitudes 

Let us start by summarizing the regularization and removal of infrared divergences. 
These arise in perturbation theory through diagrams in which virtual photons connect ex- 



ternal charged legs. The classical works of Refs. [13|,nj show how to sum the infrared 



singularities to all orders in perturbation theory and isolate an infrared finite amplitude. 
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We begin by reviewing the content of these works. Let us introduce an IR regulator A. For 
our calculation, this takes the form of a photon squared- mass, A = mt. Let A be the ampli- 
tude for a generic process involving charged particles. To all orders in the EM interaction 
A is given by the expansion 



oo 



A = J2 A * > (5) 

k=0 

with Ak = 0(a k ). Order by order one has the sequence of relations, 

Ai = a ■ aB(X) + di , 
(aB(X)) 2 

A 2 = a - ± + ai ■ aB{\) + a 2 , (6) 



Here B(\) is an infrared divergent function of A, while the are infrared finite. Summing 
to all orders results in the compact relation 

oo 

A(a) = e aB{x) -J2 a k= e aBW A(a) , (7) 

where all IR-singular dependence appears solely within the exponentiated factor aB(X), 
which multiplies the infrared finite amplitude A(a). The function aB(X) only depends on 
the external states and knows nothing on the nature of the interaction generating the process. 
On the other hand, the amplitude A(a) contains the structure dependent EM effects. These 
quantities arise naturally in the calculation described in Ref. 0, where we provided explicit 
expressions in the case of A+-. 

The construction just described needs to be supplemented by the the following comments. 
As Eq. (H) shows, the function B(X) arises as a first order correction in a. This means that 
a one-loop calculation allows resummation of the IR singularity to all orders. However, the 
definition of the IR finite part of 5(A) is not unique. This means that there is an ambiguity 
in the way one separates the infrared multiplicative factor from the structure dependent 
effects. This is a peculiar property of EM radiative corrections and does not affect the 
definition of physical observables. For example, once one picks a definition for -B(A) and 
follows it throughout the calculation, comparison with experiment will lead to unambiguous 
extraction of the EM free quantities (like ao). We will explicitly display our formulas for 
B(X) below. 

The infrared divergences of the amplitudes have now been isolated in an overall factor. 
Removal of infrared divergences from the expression for the decay rate or cross sections 
is achieved by taking into account the effect of soft real photons in the external states. 
This is motivated by the observation that for soft photons, whose energy is below some 
experimental resolution u, the generic states n and n + kj cannot be distinguished. The 
physical observable always involves an inclusive sum over the n and n + k 7 final states. We 
shall give details of this in Sect. 0. 
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B. Isospin Amplitudes 



Having described the construction of IR finite amplitudes, we can now analyze the effect 
produced on the isospin amplitudes Aq, A 2 . We start from the IR finite amplitudes in the 
charged basis {A+-, Aoo, A + o}- If is then possible to define the would-be isospin amplitudes 
by taking the following linear combinations: 

-r _ 2 1 

Aq — 2^+- + 3 00 ' 

A 2 = ^(A + --A o) , (8) 
A 2 = g^+o • 

The content of Eq. (§) is that in the absence of EM and for m u = md the amplitudes Aj truly 
describe transitions to tttc states with definite isospin. In the presence of isospin breaking 
interactions, however, the amplitudes Ai become 

~Ai = (Aj + 5Aj) e i( - Sl+ ^ , (9) 

with shifts 5A I and 7/ corresponding respectively to the modulus and phase of the original 
isospin amplitude. 

We can summarize the discussion presented so far by displaying the parameterization of 
the IR finite amplitudes in the charged basis 

A+- = (A + 5A ) e^ ^ + --*= (A 2 + 5A 2 ) e^ 2+ ^ , 

v2 

^00 = (A) + SA ) e^ 0+ ^ -V2(A 2 + 5A 2 ) e^ 2+ ^ , (10) 
A +0 = \{A 2 + 5 At) e^t) , 

where 5Aj and 7/ contain IR finite EM effects as well as other isospin breaking terms. This 
parameterization has to be compared with the isospin invariant expressions in Eq. (0). We 
recall here that this parameterization holds for the IR finite amplitude as defined in Eq. (^). 
We also observe that the shifts 5A^ and 7^ in A+q are distinct from the corresponding shifts 
in ^4 + _ and Aoo, as a consequence of the AI = 5/2 amplitude induced by electromagnetism. 
In the notation of Eq. (U) one has: 

5A = 5A 1/2 , 

5A 2 = 5A 3/2 + A 5/2 , (11) 
5 A+ = 5A 3/2 - 2/3 A 5/2 . 



C. Discussion 



The parameterization given in Eqs. (|I0|) , (|TT|) provides the basis for any phenomenological 
analysis of K — > nn decays with inclusion of isospin breaking (due to strong and electromag- 
netic effects). Comparison with experimental branching ratios (see next Section for some 
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caveats related to radiative corrections) allows one to arrive at relations among the different 
parameters. Examples of such analysis are given in Ref. [ 10 1 and the third paper of Ref. ||. 



It is legitimate at this point to ask what do we know about the parameters entering 
Eq. ( |10D and what do we want to extract from the comparison with experiment. 



1. First, let us consider the amplitudes A and A 2 . Many theoretical efforts have been 
devoted to their calculation, and no calculation can be considered fully satisfactory at 
present. We want to extract these parameters from the comparison with the branching 
ratios, eliminating the electromagnetic isospin breaking contaminations. 

2. Next, one has the shifts SA , SA 2 , and SA 2 . We have calculated the EM contributions 
to them and we are quite confident that our results capture the true underlying 
physics within the theoretical uncertainty quoted. In particular, our results provide 
an estimate for A 5 / 2 . Explicitly we find ||: 

5Aq M = (0.0253 ± 0.0072) ■ 10~ 7 M K o , 
5A™ = (0.0147 ± 0.0063) • 10~ 7 M K o , 
SA% EM = (0.008 ± 0.0088) • 10~ 7 M^o , 
A™ = (0.0137 ±0.0097) ■ 10" 7 M K o . 

For recent estimates of the size of non-EM isospin breaking effects, see Ref. ||. The 
corrections SA^f 13 ^ contain only a negligible AI = 5/2 component, as the only source 
for it would be a combination of the AI = 1 interaction proportional to (m u — m^) 
with the suppressed AI = 3/2 weak interaction. This ensures that SA^J^ft can be 
reabsorbed into A 2 . Therefore, in what follows A and A 2 still contain strong isospin 
breaking contaminations. These can be subtracted by using the results of Refs. || . 

3. Finally, let us consider the rescattering phases. In absence of isospin-breaking the 
phases 7/ vanish as a consequence of Watson's theorem. The strong phases <5/ at 
s = M\ are known through dispersive treatment of inr scattering data |. The inclusion 



of strong isospin breaking effects still gives 7/ = 0, to first order in m u — rrid |]I0" |. It 
is the inclusion of EM corrections that generates nonzero 7/. Since the phase 5 2 + 7^" 
does not enter any physical relation, we disregard it from now on. As for the phases 
7o and 72, we can relate them to EM effects in the final state interaction. This can 
be done in perturbation theory [Q] or in a more general setting provided by unitarity. 
As we shall discuss in Sect. |V|, our present knowledge of 70,2 reveals a large value 
of 72 and relies on a lowest order analysis of EM corrections to tttt scattering. Next 
to leading order corrections can be quite large at s = M^, and this adds substantial 
uncertainty to 72 . 



^For a list of references see |10| |. 
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III. ANALYSIS IN THE PRESENCE OF ELECTROMAGNETISM 



The CP-conserving sector of K — > txtx phenomenology relies essentially on three experi- 
mental numbers. These are the partial decay widths of Kg into 7r + 7r~, 7r°7r° and of K ± into 
7r ± 7T°. Knowledge of these numbers allows one to extract the invariant decay amplitudes 
and to compare with theoretical calculations. In this section we describe the procedure to 
be followed in order to extract the invariant amplitudes in presence of radiative corrections. 



A. The Fitting Procedure in Presence of Radiative Corrections 



In presence of radiative corrections, the appropriate expression to be used when compar- 
ing theory and experiment is: 



<j> " \A n \ 2 G n {uj) 



(12) 



where n = {H — , 00, +0}. In Eq. (|T2"D the left hand side Y n represents the measured partial 
width and we indicate with the parameter uj its dependence on the way soft photons are 
treated in the data analysis. It is clear that the use of different cuts leads to different values 
for the decay widths, because of the inclusion of different portions of the corresponding 
radiative channel (n + 7) in the data sample. On the right hand side of Eq. ( |P2"D one 
has the kinematical parameters s n (squared center-of-mass energy of the process) and <3> n 
(the two body invariant phase space associated with the decay). The quantities related to 
the dynamics are A n and G n (u). A n is the infrared finite invariant amplitude, as defined in 
Sect. [II], Eq. fllPp. It contains the true weak transition component that we wish to ultimately 
extract as well as infrared finite electromagnetic corrections. G n (uj) is the infrared factor 
associated with the combined effect of virtual and real photons. The latter contribution 
involves an integration over the soft-photon phase space: this has to be done with the same 
prescription used in extracting the experimental number T n (u). 

The extraction of IR-free invariant amplitudes is a straightforward consequence of 
Eq. (12). Specializing to the K — ► Tin case, one has: 



\A, 



M K o 

\Aqo\ 

M K o 
\A +0 \ 



where 



?7T T + _(CJ) 



y/2 \ P+- 




(13) 



M K + \ p +0 G +0 (uj) 



P+- 



M 2 4 
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Poo 



P+o 



( M K , 



MK 



(14) 



+ 



M 2 - M 2 + ' 

7T U 7T + 

2M K + 



-M_ 2 „ 



In order to carry out the procedure one needs the physical masses, the experimental input 
r n (c<j), and the corresponding theoretical infrared parameters G + _ )+0 (co>). If things are done 
properly, the u dependence cancels in the ratios on the right hand side of Eq. which 
provides then the values of \A n \- Having the amplitudes \A„.\ one can then use Eq. ([U]) to 
obtain relations between the physical parameters. However, there is an important issue to 
be addressed in order to complete the program outlined above: it is the proper definition of 
the infrared factors G H (u>) and G + o(uj), to which we now turn. 



B. The Infrared Factors 

As already noted, the theoretical definition of G + -(u) and G + o(u) involves integration 
over the soft-photon phase space, with cuts generically indicated by to. The infrared factor 
is experiment dependent and accounts for the component of the radiative mode (irn'-f in our 
case) included in the parent mode (7T7t) branching ratio. We shall work to order a and thus 
include only the effect of the incy radiative state. While we display the ingredients needed 
for any experimental analysis, we carry the treatment to completion for the case where 
the infrared sensitivity is isotropic in the center of mass. That is, we display the explicit 
formulas obtained when one integrates over the nn^ phase space applying an isotropic cutoff 
uj on the photon energy E 1 in the center of mass frame. This is most natural to use in the 
data analysis for an experiment such as KLOE at DA$NE, given the working conditions 
of the machine and the detector geometry. A more detailed study (either theoretical or 
experimental) is required in order to apply radiative corrections to the data analysis of other 
experiments. The presence of several high statistics experiments offers a unique opportunity 
for performing an accurate measurement of K — > tttt branching ratios, including the effect 
of radiative corrections. Such an analysis would fill the present gap in the study of radiative 
corrections to K° — > tttt, and would therefore be highly desirable. 

Calculation of the factors G + ^(uj) and G + q(uj) requires consideration of a combination 
of effects due to virtual photons in the amplitudes A+- i+ o (field theoretic version of the 
non-relativistic Gamow factors) and to soft real photons entering the process K — > tctt + 7. 
Moreover, in the infrared region of the spectrum, the amplitude for K — > 7T7T+7 is dominated 
by the internal bremsstrahlung (IB) component, proportional to the non-radiative amplitude. 
For the case in question one has: 



A IB 


= eA+- 




■p+ 


e 


■P-\ 


! 


■p+ 


Q 


■P-) 


A IB 


= eA +0 | 




p+ 


e ■ 


Pk\ 




p+ 


Q 


Pk) 



(15) 



where e and q are the polarization and momentum of the emitted photon. Now, the infrared- 
finite observable decay rate is 
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[00 



(16) 



where n = H — , +0 and 



vw . d$ n \A n (l + a5 n (m 7 )) | , 
zivix 



ZMk JE~,<u) 



2M, 



K 



d$ n \A n \ 2 I n (m 7 ,a/) 



(17) 
(18) 



At this stage the rest of the calculation becomes dependent on the geometry of the exper- 
iment. For a kaon at rest with an isotropic detector, the acceptance cutoff E 1 < uj will be 
the same in all directions. For a kaon in flight, the acceptance may be different for photons 
emitted in different directions. In this latter situation, the integral in Eq. ( ]T8[ ) needs to be 
numerically integrated over the detector acceptance and then added to the two body result, 
Eq. (|17])- This will then be finite in the limit of m 7 — > 0, and will allow the measurement of 
the IR finite amplitude A n . We carry out this procedure explicitly below for the case of an 
isotropic cutoff. 

In Eqs. (|IT|),(fi"8|) d§ k is the differential phase space factor for each process, A n is the 
IR finite amplitude, as defined in Sect. || by extracting the IR divergent functions B n . The 
B n can be calculated by considering one-loop diagrams with virtual photons connecting the 
external legs and using a point-like vertex for the weak interaction. The definition of B n 
is not unique due to the possibility of adding IR and UV finite terms to it. The explicit 
expressions given below | fully specify our choice of what goes into B n and what goes into the 
structure dependent shifts 5A n . The second expression in Eq. (|18D describes the factorization 
of the Ten and 7 phase spaces, valid with an accuracy of uj/Mk- Explicitly one has: 



d 3 q 



e^u, (2tt) 3 2E 1 



E 

pol 



AIB 2 



A 



Combining the above results one arrives at 



G n (uo) = 1 + 2aTleB n (m J ) + I n {m 1 ,u) + 0(a 2 ) . 

We now collect the explicit form of the functions _B + _ j+ o and 1^ 
definition of G + _ ;+ o (see Eq. |2(|). For G + -(u) one needs: 



(19) 
(20) 

+ o, entering in the 



J + _(m 7 , lo) 



1 

47T 

O 

7T 



mi 



ft 



mi 



■w>-(£) ! 



(21) 



where 



' See also Ref. 
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P = (1 -AM 2 / 'Ml) 1 ' 2 



(22) 



and 



a(/3 ) = l + l±^ln^-^ 



1 + /5 2 

2(3 



2(3 \l + (3) 

1 + /5, 1-/5 2 



it 2 + In 



In 



1 -(3 4/5 



+ 2/ 



1 + 

2/5 , 



2/ 



, 2/5 , 



+ 2 + /51n 



l + (3 
l-(3 ' 



1 1 + /5 1 + /5 
— In H — 



/5 1-/9 2/5 



2/(-/5) - 2/(0) + / 



'1 + /T 



-/ 



'1-/T 



2 1-/5 



ln(l -/5 2 ) + ln2 In 



1-/3 
1 + /5 



/(*) = -f 
Jo 



x ~y 

dt — In II — tj 

t 1 1 



(23) 



We note that Z? H includes the Coulomb factor -Kct/v TC \ (first term in H + _((3)), as well as 

typical field theoretic effects. As for G +0 (ui), one has: 



B +0 (m 2 ) 



I + o(m-y,uj) 



1 

47T 

a 

7T 



2b{(3) \ n ^f + H +0 ((3) 



H/3)ln(^) 2 + F +0 (/5) 



(24) 



where 



and 



1 (l-(3 s 



#+o(/?) = 



1 

7 



log 



1 + /?\ 



log 



-2 log 



'l-(3' 



+ log 



+ / 



4/5 \ 



log 



log 



2 + 2/5 s 
2-2/5 \ 



'l+/5^ 



+ / 



+2 + log 



4/5 

1-/3 2 



-/ 



1-/5 

4/5 



1 + /5 N 

. 2/5 / 



/ 



7?-r 

, 2/5 , 



2 ^ ^ 



1-/5 

'1 + /3 S 
.1-^ 



2 ^ ^ 



1 + /5 



11 



D(x) 

E(x) 
p(x) 



l-(3 2 

(x — Xi)(x — X2 



dx — — log 

-l D(x)p(x) 



E(x) + p(x) 
E(x) -p{x) t 



Ml 
Ml 
Mk 

4 
Mk 

4 



1±(3)-1 

(3-x) 
P(l + x) . 



(25) 



G + _ and G +0 do not depend on the infrared regulator m 7 . We display plots of the 
functions G + ^(u) and G +0 (u) on a typical range of values for the parameter u in Figs. || and 
|3|. The functions G n (u) are the only EM effects previously considered in the literature j5|, 
although with a slightly different definition. In fact, the works of Ref. H use a point- 
like vertex for the weak interaction, and therefore are not sensitive to structure dependent 
corrections. In these works the EM effects due to wavefunction renormalization and vertex 
correction go entirely in the definition of B n (this includes also the UV divergent terms, 
regulated by means of a cutoff). Apart from the cutoff-dependent term and an extra finite 
contribution, our expressions match the ones given in the second and third papers of Ref. || . 
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FIG. 2. The function 



IV. EFFECT OF EM ON K -» vrvr PHASES 

The amplitude parameterization we have used for the previous analysis already implies 
that K — ► 7T7T data do not provide direct information on the strong tttt phase shift difference 
8q — 62- This would only be true in the isospin limit (77- = 0). In this limit, the requirement 
that strong interaction phases appear in the weak decays is known as Watson's theorem, 
which is valid whenever the final state rescattering involves only elastic scattering. Despite 
the fact that isospin breaking has been long understood to cause mixing of weak amplitudes, 
there has been no recognition that the strong interaction phases no longer suffice to describe 
the rescattering effects. This occurs because elastic rescattering is no longer the full content 



12 



G+o(") 



1 .01 



w (MeV) 



5 



10 



15 



20 



. 99 




. 98 



. 97 



. 96 



FIG. 3. The function G+o(u). 



of final state interaction, so that the conditions for the application of Watson's theorem no 
longer apply. Moreover, there exists a sizeable discrepancy between the determination of 
5q — 62 from K — > ttti data using isospin relations and the favored value of the phase shift 
difference known by other determinations. This seems to point to violations of Watson's 
theorem. It is our purpose to set the framework for the correct treatment of this problem 
in the isospin breaking real world. We shall accomplish this by writing a coupled channel 
unitarity constraint in the presence of EM interactions (and isospin breaking in general) 
and solving for the parameters 70 and 72 entering Eq. fllDp. This analysis will complement 
and extend the perturbative results obtained in Ref. 0, which already indicated a large 
value of 72. We defer to the next section the extraction of So — 62 and the uncertainty to be 
associated with it. 



The first step in our program is writing down meaningful unitarity relations in the 
presence of EM interactions. Here, it is natural to work in the charged basis {tc + -k~, 7t°7T } 
and then try to recover the notion of isospin amplitudes. In order to fix the notation, let us 
start from the unitarity relations involving the decay amplitudes of K° to {tt + tt~, 7r°7r } in 
the limit in which EM is turned off. Then, only the tttt intermediate states have to be taken 
into account and one finds: 



In Eq. A^ and ^4oo represent the K° decay amplitudes and 7j ;i is the T-matrix element 

for the transition % — > / ( in this case it only involves pion-pion scattering ). 'T* x A' denotes 
the product of amplitudes integrated over the intermediate state phase space. In the case 
considered here of two-pion intermediate states, one has: 



A. Extended Unitarity Relations 




% (t;_.+_ x A + - + T * 0;+ _ x Aoo) 




(26) 




(27) 
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where (3 = (1 — 4m 2 r /m| : )™ 1 / 2 is the pion velocity in the kaon rest frame. $ s is the symmetry 
factor for identical particles (equal to 1/2 for the 7r°7r° state) and T is the S-wave projection 
of the TTii scattering amplitude T(cos#), defined by: 

T = f +1 d(cos9)T(cos9) . (28) 

647T J -I 

Turning on EM interactions introduces isospin breaking dynamics as well as IR singulari- 
ties in the amplitudes and the opening of intermediate radiative channels. Specifically, A+-, 
7oo,+-, and X|__ >+ _ become IR divergent and 7^_ j+ _ acquires a purely Coulomb component 
(also IR singular). The work of Refs. jl3lJl4jl , summarized in Sect. [II] teaches us that one can 
always isolate the singularity in a multiplicative exponential factor 

Af ti = e aB f- A u . (29) 

Here Af,i is the IR finite amplitude and Bf^im^) is the IR singular factor that depends 

only on the external states. We shall only need the factor , already encountered in this 

paper, associated with a pair of charged pions in the initial or final state. We note that 
5j j(m 7 ) is in general complex. In particular, its imaginary part is equal to the Coulomb 
scattering phase shift associated with each pair of charged particles in the initial and final 
states |l4f . 



Upon integrating over the phase space and using the above mentioned property on the 
Coulomb phases, one can rewrite Eq. (|26|) in terms of IR finite quantities and Bfj factors. 
Moreover, due to the integration over the phase space some contributions in the Bfj factors 
simplify and one ends up with: 



A+- -X + _=i (t;_ ;+ _ x A+- e 2aV - eB+ ~ + r m+ _ x Ac, 
Aoo - A* 00 = i (r;_. 00 x A + - e 2 ^~ + T 00]00 x A 00 ) . (30) 



Note that now T + _ + _ is the IR finite 7r + 7r~ — > 7i + ti~ amplitude subtracted of its purely 
Coulomb term. 

Eq. ( PD| ) contains IR singularities, but the analysis is still missing an important effect 
of EM: the opening of inelastic radiative channels. This is the key ingredient in obtaining 
an IR finite set of unitarity constraints, as it was in obtaining an IR finite cross section or 
decay rate. In fact, the IR singularities will cancel in the sum over the n + ir~ and 7r + 7r~7 
intermediate states, with the same mechanism described in the definition of r +0 earlier 
in Sect. PTT| Working at order 0(a), we consider only the radiative state 7r + 7r - 7. For our 
analysis we require the amplitudes for K° — > 7r + 7r~7 and tttt — > 7r + 7r~7. We include only 
the internal bremsstrahlung component of these amplitudes, known to be dominant over 
possible direct emission terms. Now one has to integrate over the full 7r + 7r~7 phase space 
and the final result for the unitarity condition reads: 

Im So ) = K 2 t + : ; V(i 1 + + A A + + _ ) ) tZo ) So ) ■ (31) 

We recall that T a j> are the S-wave projections of the tttt scattering matrix. A + _ is the IR 
finite remnant of the sum of IR singular terms in the 7r + 7r~ and 7r + 7r~7 intermediate states. 
In terms of the notation of Sect. |111 BL it is given by: 
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A, 



25MI 



+ 2aKeB + _(m 1 ) + e 2 



E 

pol 



9- 



• k 



(32) 



Here the first term is the phase space correction due to the EM mass-shift of charged 
pions. The second term is the effect of infrared virtual photons, while the third term is the 
effect of real soft photons in the intermediate state 7r + 7r~7. Numerically we find (displaying 
separately the phase space contribution and the remainder): 



A 



+- 



-14.8 + 10.8) ■ 1(T 3 = -4.0 • 10' 



(33) 



B. From Charge to Isospin Basis 

Assuming unitarity of the S matrix, we have thus far obtained a set of relations containing 
the IR finite amplitudes in the charge basis. In order to compare with usual treatments of 
this problem, we rotate now to the isospin basis for the K — > ittt amplitudes, 



3iso= 7 =; .% \m m ■ (34) 





Applying the same transformation to the whole system in Eq. Q3~ID one obtains in matrix 
form: 

XmA lso = $ (T\ so +n)A lS o , (35) 

where 

T ISO =(l° I 02 ) , (36) 



T20 T2 



and 



_ / Rqo R02 \ _ 1 A ( 2T * y/2T* \ 
^-[R^o R22 ) ~ 3 + - { V2T* T* ) ■ (37) 

Tiso is the tttt scattering matrix in the isospin basis, while the matrix TZ, proportional to 
A + _, contains the effect of IR radiative corrections and the radiative intermediate channel. 
The 7T7T scattering T-matrix now involves both strong and EM interactions, and thus contains 
isospin-violating matrix elements. In the conventions used in our work, the amplitudes for 
the 7T7T scattering in the isospin basis are expressed in terms of the charged ones as: 

Tq — - (4T + _ i+ _ + T 00] oo + 4T 00j+ 

T2 = - (T+- t +- + T o,oo — 2T o,+-) , (38) 

?20 = T 2 = (2T + _ + _ — T o,oo — ^00, 
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A general parameterization of the irn transition matrix in the isospin basis is: 

(rio e 2i5 » - l)/(2i) ae^ 5o+s ^ \ , . 

ae i(5 +S 2 +A) ( m e 2iS 2 _ l)/(2i) J • [6J) 

In this parameterization we allow for isospin mixing (the off-diagonal parameter a) and for 
possible non-unitarity in the tttt two dimensional subspace (due to opening of other channels). 
This is accomplished by introducing the inelasticity parameters ^0,2 arid the extra phase A 
in the off-diagonal term. The parameters rjj are of order 1 + 0(a 2 ), while a and A are of 
order a. 

The form given in Eq. ( |39| ) is fully general and includes all isospin breaking effects. 
However, strong isospin breaking is expected to induce only subleading rescattering effects. 
In fact, m u — m^ 7^ produces an I = 1 perturbation to the original interaction. This is not 
sufficient to mix the 1 = and I = 2 tttt scattering states when treated to first order, nor 
does it cause a splitting of the masses of the charged and neutral pions. This implies that 
elastic scattering of these states is still the only option, and to first order in m u — rrid the 
parameter a does not receive contributions. The values of the phases 5i may in principle be 
slightly modified by the quark mass effect, yet this is contained in the measured values of 
the experimental phase shifts. 




C. Solution for 70 2 

We are now in position to explore the consequences of unitarity on the rescattering phases 
7o,2- We write 

A I = A I e i{Sl+ ~ fl) (40) 

and insert these expressions into Eq. (^). We then solve for sin 70 and sin 72 to first order 
in a, taking into account the AI = 1/2 hierarchy of magnitudes. After some simple algebra, 
we obtain the solutions 



shi7o = (3 (TZeRoo — tan5 Xmi?oo) — O(asin5 ) 

(41) 



sin 72 = P=r 

A 2 



T20 H r (R-eR2o cos S - 2mR 2 o sin 5 ) 

COSO2 



The most important feature of these results is the factor A /A 2 in the formula for sin 72. 
This implies that even though the non-elastic scattering is electromagnetic in origin, it is 
enhanced by a large factor that allows the net change to be significant. Eq. ([11]) gives us the 
desired expression relating the phase 72 to isospin breaking rescattering effects. These are 
contained in the parameters T 20 , the mixing amplitude between tttt states, and R 2 o- This 
last parameter contains the effect of the radiative intermediate channel 7r + 7r~7 as well as the 
phase space correction. We note that Eq. ( f4*l"D is a generalization of the relation obtained at 
one loop in ChPT. However, inspection reveals that the perturbative determination contains 
only the phase space effect and the T 2 o mixing in lowest order. 

In attempting to estimate the magnitudes of the new phases, we are hampered by the 
fact that the analysis of electromagnetic effects in tctt scattering is not yet complete in the 
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literature. Two groups have provided analyses of reactions involving neutral mesons |16 
but the channels with all charged particles are not yet fully analyzed. We require the 
scattering elements at center-of-mass energy equal to the kaon mass. The threshold matrix 
elements are known from simple tree level calculations, and we will use these in our estimate 
below. However, the amplitudes can experience large changes at s — M\, and one needs at 
least one-loop chiral perturbation theory in order to obtain these. As these results become 
available, they can be used to update our numerical estimates. 

We estimate the off diagonal parameter at lowest order in chiral symmetry obtaining: 

T 02 = ^ . !Ml ^ 2 . 7 x 10 -3 . (42) 
For the parameter R 20 , proportional to the radiative effect, one has 

R 20 = ^A+_T 2 * , (43) 

and we use the form 

T 2 = ie* fo sin5 2 , (44) 

with the phenomenological central value of 6 2 = —7.0° [?]. Numerically this leads to 

TZeR 20 = 0.280 • 10~ 3 

ImR 20 = 0.034 ■ 10" 3 (45) 

These numerical estimates allow us to identify the off-diagonal AI = 2 rescattering as the 
major new ingredient in the final state phases and to arrive at the result: 

7o = -0.1°, 72 = 3.1° . (46) 

We note here that the result for 7 2 is quite large, amounting to almost 50% of the strong 
phase 5 2 at s = M\. 



V. SAMPLE FIT TO K -» vrvr DATA 

In this section, we provide a tentative fit to the present experimental data. This is meant 
as an illustration of the ideas that we have discussed above, and hopefully will provide a 
model for a new fully consistent experimental analysis of new data, taken with the full 
treatment of electromagnetic effects. We describe our treatment as tentative because it 
involves older data sets which were taken without the inclusion of radiative corrections. We 
cannot fully account for the experimental acceptances, and are forced to adopt a cruder 
procedure. However, the sample fit is none the less of interest because it illustrates the 
significant sensitivity of various quantities to electromagnetic corrections, and represents 
the best that can be done with the present data set. 
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A. Data Analysis 



It will be convenient in the discussion to follow to first define 

X% = $i + li {i = 0, 2) . 



(47) 



Then having G + -(tu), G + o(u) and the structure dependent corrections 5Af , one is in a 



position by using Eqs. (|13[) and (|T0| ) to extract the quantities A , A 2 and Xo — X2- 
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FIG. 4. Fitted Aq as a function of u; 
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FIG. 5. Fitted Ao as a function of u 



As experimental input for the branching ratios we use the PDG averages, although these 
numbers come with no reference to what portion of the K — > 7m / ~f mode is included. In order 
to understand the attendant uncertainties and ambiguities of this approach, in Figs. £|, || 
and H we plot the output of our fit as a function of the parameter u>, the upper cutoff 
for IR photons in the center of mass frame. It is not clear to which value of uj (if any) 
the experimental numbers correspond. This ignorance gives rise to little uncertainty in 
the extraction A 2 and to a moderate one in the extraction of Aq, for u varying between 
1 MeV and 20 MeV.Q However, more delicate is the situation for the extraction of the 



**This range is chosen to reflect a realistic possibility for detector resolution. 
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FIG. 6. Fitted xo ~~ X2 as a function of co . 



phase xo — X2, where a variation of the order of 10% is seen over the considered range of 
u. Thus our analysis indicates that the extraction of rescattering phases from K — > 
data is sensitive to the treatment of soft photons. In the absence of precise experimental 
information, it is not possible to pick a definite central value for our output. We thus quote 
the results for the set of EM-free quantities with two error bars. The first one is due to 
the spread in the central values according to variations of u between 1 and 20 MeV. The 
second one comes from propagating the experimental uncertainty in the decay widths and 
the theoretical uncertainty on the inputs 5Aj M . We find: 



Xo 



A = (5.450 ± 0.020 ± 0.015) x 10" 7 M K o 
A 2 = (0.255 ± 0.001 ± 0.009) x lQ~ 7 M K o 
■ X 2 = (56 ± 4 ± 4)° . 



(48) 



These results should be compared with the ones presented in Eq. ([|), derived from the 
analysis in the isospin limit. The most important new feature is that considering EM 
corrections places larger error bars on all these quantities. In the case of A2 the reason 
for this resides in the quite large theoretical uncertainty on SA^. For A and the phase 
difference Xo~ X2, the larger error bar is due essentially to incomplete information concerning 
the treatment of the radiative channel. A measurement of the the partial width r + _(co>), 
with accuracy level of ~ 0.5% (this is the accuracy level of the present PDG numbers), 
accompanied by information on soft-photon cuts, would allow one to extract a definite 
central value for A and Xo ~ X2- As a consequence, this would eliminate the first error bar 
associated with A and xo ~ X2 in Eq- (|48|), reducing the total uncertainty by 50% or more. 



Indeed, such an analysis will be performed by the KLOE experiment at Da$ne |L5 . 



B. Extraction of 5o — 62- Discussion 
Finally we turn to the extraction of 5q — 62 from K —>■ tttt data. The relation to be used 



is: 



S0-S2 = (xo - X2) &t + 72 - 7o • 



(49) 



19 



As shown by Eq. (f49[), the extraction of the strong phase difference relies on two distinct 
inputs: 

f . The first one comes from the fit to K — > ittt branching ratios, which provides Xo ~ Xi- 
In Sect. |V|A we discussed such a fit and pointed out the sensitivity of xo ~ X2 to 
cuts used for soft real photons. In the absence of information on these cuts a precise 
determination of xo ~ X2 is not possible and our conclusion is that the error bars are 
larger than previously thought: 

Xo ~ X2 = (56 ± 8)° (50) 



2. The second input concerns the magnitude of 72 — 70- We have thus far established 
a general framework (based on unitarity) for the analysis of these isospin breaking 
phases. We found that 72 receives a A J = 1/2 enhancement and the dominant effect 
is due to isospin mixing in the irn rescattering rather than to radiative intermediate 
channels. We have provided an estimate of T 2 o at lowest order in the chiral expansion, 
leading us to write: 

72 - 7o = 3.2° + 7 f p2) • (51) 

(e 2 v 2 ) 

The possibility of large chiral corrections to T 2 o (associated with 72 ) cannot be ruled 
out, given the results obtained in the analysis of other EM corrections (violations of 
Dashen's theorem and the K —>■ tttt amplitudes). 

In light of the previous discussion, we give the following value for 5 — 5 2 from K — > nn data: 

5o-5 2 = (59 + 7? V) ±8)° . (52) 

The leading order estimate for 72 is seen to worsen the discrepancy between the central values 
of weak and strong determinations of Sq — 5 2 - However, the large uncertainty associated with 
radiative corrections makes impossible a precise comparison at this stage. In this sense, the 
phase puzzle is alleviated, its cause being a previous underestimate of error bars. Indeed we 

(e 2 p 2 ) 

believe that the combined effect of radiative corrections to xo — X2 and calculation of 72 
can fully resolve the puzzle, providing a satisfactory theoretical formulation of the problem. 
In fact, once a more precise extraction of Xo — X2 becomes available, Eq. can be used 
to extract T20, and thus information on the isospin breaking dynamics in irn scattering at 



VI. IMPACT ON CP PHENOMENOLOGY 

In the present section we focus on the consequences of our work to CP phenomenology in 
the kaon system. Our work gives rise to interesting effects only in the theoretical analysis of 
e'. In particular, we provide an estimate of the isospin breaking parameter f2 EM , the effect 
of the A J = 5/2 amplitude and the phase of e'. 
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The analysis of direct CP-violation in K — > txtx proceeds exactly as in the standard 
case, except that now we work with the IR finite isospin amplitudes Aj and the final state 
interaction phases xi associated with them. One can then write 



y/2 



J(X2~Xo) 



TleA 2 
TleAr 



ImAo I171A2 



TleAr 



TZeA, 



Defining 



TleA 2 



UJ 



TleA 

and neglecting the small effect of SA /A one arrives at 

XtuAq 



e = — 



V2 



e i(X2-Xo) u 



TleA 



1 ZmA 2 
1 - 

uj XijiAq 



(53) 



(54) 



(55) 



We recall here that in the Standard Model analysis the imaginary part of A is generated by 
the so called gluonic penguin, while the phase of A 2 is generated by the electroweak penguin. 

In order to make manifest the effects of electromagnetic corrections, we now further study 
Eq. ( |55|) . The first new effect is to be found in the parameter To. It is due to the presence of 
the AI = 5/2 amplitude, distinguishing A 2 from A\ (see Eq. fllDD). In the usual treatment 
one uses the parameter 



TleA 



to 



TleAo 22.2 

However, our derivation shows that one should use uj. The two are related by: 

_ TleAz TleA 2 



uj 



TleA JZeA 



=+=w(l + fs/2 



(56) 



(57) 



The other relevant phenomenon is the leakage of the octet amplitude into A 2 , providing 
the dominant part of 5A 2 . This brings an extra contribution to the CP- violating phase of 
A 2 , essentially generated by the gluonic penguin and transferred to A 2 via isospin breaking 
effects. This mechanism is usually parameterized by: 



n 



iso— brk 



1 lm5Al?°- bTk 



to 



XmA 

where f2 lso ~ brk will have contributions from both electromagnetic effects (f2 EM ) 
strong interaction effects (f2 STR ) associated with m u 7^ m^, 

fi iso-brk = ^EM + fi STR _ 



The above observations lead us to write: 
V2 TleA 



e = — - 



1 - 



1 lmA 2 
uj ImA 



+ fs/2 - n 



iso— brk 



(58) 
and from 

(59) 
(60) 



Comparing Eq. (|60|) to the standard analysis (not including EM corrections), one identifies 
three new effects. 
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1. The factor / 5 / 2 appears: it can be obtained by inserting in Eq. (|57|) our previous 
estimates of SA 2 and SA 2 (see Ref. [|]). We find 

/™ = (9.3±6.1)-l(r 2 . (61) 

The large uncertainty reflects the one in SA 2 . We thus find that this effect tends to 
increase (although slightly) the central value of e'/e. The authors of Ref. |K| find an 



opposite result because they use the "phenomenological" value of A 5 / 2 . We believe 
that the phenomenological determination of A 5 / 2 , as performed in Ref. ||10|| , suffers 



from large systematic uncertainties due to neglecting IR effects and the EM phases 7/. 

2. One has to consider the electromagnetic contribution f2 EM , to be added to existing 
estimates of f2 STR due to strong isospin breaking. Again, the analysis performed in 
Ref. [|J enables us to get the magnitude of Q EM , since we calculated there the octet 
induced component of 5A EM . Thus we can write: 

^ EM _ Tie A _ lmSA 2 _ Tie A _ Tie SA 2 

~ TieA 2 ' lmA ~ TleA 2 ' TieA ' ^ ' 

Numerically we find: 

fi EM = (6.0 ±2.5) ■ 1(T 2 . (63) 



3. One observes that the phase of e'/e is related to Xo ~ Xi an d not to 5 — 5 2 , although 
with the present accuracy it is hard to make a meaningful determination. We find 



^' A =(x 2 -Xo + 0-^ = -(ll±8)° . (64) 

The resulting effect on the real and imaginary part of e'/e is below the sensitivity of 
present kaon factories. 

We conclude by observing that the individual terms f§/ 2 and Q EM have a respectable size 
but enter in the expression for e' with opposite sign. The net effect has a very small central 
value with a large uncertainty. 



VII. CONCLUSIONS 

In this paper we have attempted a full phenomenological analysis of K — > tttt decays 
in the presence of electromagnetic interactions. We have provided a general parameteriza- 
tion of K — » 7T7T amplitudes to include the effect of isospin breaking interactions. Such a 
parameterization has allowed us to organize the calculation in terms of three main effects: 
structure dependent corrections (see Refs. fl|-|3|]), electromagnetic infrared corrections, and 
isospin breaking in final state interactions. We have also studied the effect of electromagnetic 
corrections on the direct CP- violation parameter e'. 
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A. IR Effects: Need for New B.R. Measurements 

It is well known that the calculation of IR effects requires knowledge of the experimen- 
tal cuts used in treating the soft photons emitted in the K —> mr decays. In Sect. [V] we 
have pointed out that the PDG numbers come with no information concerning the radia- 
tive channel, and this seriously compromises any attempt to properly include the radiative 
corrections. In the absence of experimental input, we have performed a calculation of the 
IR effects in a simple theoretical scheme (isotropic cut on the photon energy in the center 
of mass system). We have shown how this incomplete state of affairs produces uncertainties 
larger than previously thought in the EM-free quantities 

We strongly urge that a measurement of the K — > branching ratios be performed at 
one of the current high statistics kaon experiments. To be precise, it would be interesting to 
have a set of measurements of T n (u) (n = H — , +0) at different values of u (the soft photon 
upper cutoff in the center of mass frame). This would allow anyone to apply our calculation 
of G + _ and G +0 in making a phenomenological analysis (as in Sect. |V|). Of course, each 
distinct experimental procedure would require its own theoretical calculation of G + _ i+0 . All 
such studies would be equally welcome, as long as they provide information on the inclusive 
sum of TT7T and tttcj channels. We stress that such measurements are necessary in order to 
fully address the impact of EM on K — > ttti decays. 

B. Final State Interaction Phases 

We have shown that isospin breaking changes the description of rescattering phases in 
K — > 7T7T decays, as Watson's theorem is no longer applicable. We have described this new 
feature within the general framework provided by the unitarity relations, pointing out that 
the relevant effect is of electromagnetic origin. In Sect. [TV] we have set up the framework 
relating the extra phases 70,2 to EM effects in nir scattering. Our leading order analysis 
finds a large effect in 72, equal to 50% of the strong phase 82 ■ The general framework 
presented has the potential to fully resolve the long standing inconsistency between the 
strong determination of 5q — 82 at s = M\ and the one emerging from K — > ttti data. At 
present, little can be concluded due to the large uncertainty in the phase xo ~ X2 and the 

(e 2 p 2 ) 

lack of a calculation for 72 . The first problem will be solved by new measurements of the 
branching ratios (including proper information on radiative effects). The second problem 
depends on the theoretical ability to calculate EM corrections to tctt scattering at order e 2 p 2 
in the chiral expansion. 

C. CP Phenomenology 

Finally, we have analyzed the impact of electromagnetic corrections on CP phenomenol- 
ogy (see Sect. |VJ), pointing out the new features in the study of e'/e. The isospin breaking 
effects can be encoded into the factors Q and f^/2, and also affect the phase of e'. Both f 5 / 2 
and Q receive contributions from strong isospin breaking and electromagnetism. We have 
provided an estimate for the electromagnetic effect, finding results of the order of 10%, for 
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these parameters. They appear with opposite sign, and thus do not produce sizeable shifts 
in the theoretical prediction of e'/e. 
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